Hole pairing and phonon dynamics 
in generalized 2D t — J— Holstein models 
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The formation of hole pairs in the planar t-J model is studied in the presence of independent 
dynamic vibrations of the in-plane oxygen atoms. In-plane (breathing modes) and out-of-plane 
(buckling modes) displacements are considered. We find strong evidences in favor of a stabilization 
of the two hole bound pair by out-of-plane vibrations of the in-plane oxygens. On the contrary, 
the breathing modes weaken the binding energy of the hole pair. These results are discussed in the 
context of the superconducting cuprates. 



PACS Numbers: 71.27.+&, 71.38.+i, 74.20.Mn 



The electron-phonon interaction plays the key role in 
the conventional BCS theory of superconductivity. It 
is the source of the effective (retarded) attraction be- 
tween the electrons and hence of the dynamical effect 
for pair formation. On the contrary, in unconventional 
superconductors like the high-Tc cuprates, the driving 
force for superconductivity is commonly believed to be 
the strong electronic correlations. However, it is theoret- 
ically known that, in strongly correlated systems, even 
moderate electron-phonon interactions can have drastic 
consequences. For example, it can enhance charge den- 
sity wave (CDW) and spin density wave (SDW) insta- 
bilities due to polaronic self-localization effect Ex- 
perimentally, the observation of some oxygen isotope ef- 
fect in the high-Tc cuprates [|j has given evidences for 
some contribution of the electron-phonon interaction in 
the superconductivity, even though the dominant pair- 
ing mechanism is due to strong antiferromagnetic corre- 
lations. The interplay between strong electronic corre- 
lations and electron-phonon interaction still remains an 
open question. 
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FIG. 1. Schematic lattice displacements of breathing (a) 
and buckling (b) modes in the Cu02 plane. 

For the sake of simplicity, we describe here the low 
energy electronic degrees of freedom by a single band 



t-J model. We also restrict ourselves to the vibrations 
of the in-plane oxygen atoms of the Cu02 plane which 
have been shown to be essential. Two types of displace- 
ment have to be considered: (a) in-plane breathing modes 
and (b) bukling modes, as shown schematically in Fig. 
|l|. When the equilibrium position of the oxygen atom 
lies away from the Cu plane by uq in Fig. |l|(b), the 
electron-phonon interaction becomes linear in the oxy- 
gen displacement perpendicular to the plane Q], as it is 
always the case for the breathing modes. Such a buck- 
ling structure is realized in YBa2Cu307_i-. In the anti- 
adiabatic limit the two modes of interaction give an ef- 
fective nearest-neighbor (NN) hole-hole repulsion (a) and 
attraction (b) respectively. In a week coupling i-matrix 
approximation which included an RPA antiferromagnon 
spin-fluctuation exchange and a phonon exchange, Bu- 
lut and Scalapino found that the buckling mode 
can enhance dx2^y2 pairing. Using an antiferromag- 
netic induced hole dispersion and treating the electron- 
phonon interaction at the mean-field level, Nazarenko 
and Dagotto found that the buckling mode can give 
rise to a (1^2 _y2 wave superconducting ground state (GS). 
However, both of these results involve uncontrolled ap- 
proximations which are inadequate for treating the Hub- 
bard and t-J models in the absence of phonons. Thus 
it is of interest to carry out a numerical investigation of 
this problem. Our results are based on exact diagonaliza- 
tion studies of small t-J-phonon clusters. In agreement 
with the approximate results we find that the breath- 
ing mode suppresses the two- hole pairing while the 
buckling mode stabilizes it |^,^. However, in addition, 
we have examined the effect of the phonons on the kinetic 
energy, antiferromagnetic structure factor and hole-hole 
correlations, giving a more detailed picture of the role of 
dynamic lattice vibrations on the hole pairing. 

The hamiltonian is a generalization of the t-J-Holstein 
Hamiltonian M, 
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where is the usual hole creation operator, ni and 

are the electron and hole local densities respectively, 
m is the oxygen ion mass, £7 is the phonon frequency 
and (5 = X, y differentiates the bonds along the x- and 
y-direction respectively. The sign - (+) in the last term 
corresponds to the breathing (buckling) mode. Through- 
out, energies are measured in unit of the hopping inte- 
gral t. The electron- phonon g-term involves the coupling 
of each copper hole with the displacements of the four 
neighboring oxygens ui.s and This is clearly dif- 

ferent from the on-site Holstein coupling ^ which has 
been recently used to mimic the coupling with the api- 
cal oxygen modes in the framework of the t-J model . 
Note that the displacements ui,s a-re considered through- 
out as independent variables. 



lattice with Uph up to 5 |lO|] . The behaviors with the cou- 
pling constant Aq of the various relevant physical quanti- 
ties are found to be unsensitive to Uph- However, Uph = 1 
generally underestimates the role of the phonons. Clearly 
the one-phonon calculation is a good approximation in 
the weak-coupling region (J = 0.3,0 = 0.2, A < 0.3 ). 
This truncation procedure enables us to study \/8 x a/8 
cluster with all the phonon modes (16 modes). We inves- 
tigate the one and two-hole GS of hamiltonian (|l|) in a 
regime (0.3< J < 0.5) where, in the absence of phonons, 
the two-hole pairing state is stabilized by the antiferro- 
magnetic correlation, and we take a realistic phonon fre- 
quency = 0.2. Since the a/S x a/S cluster with periodic 
boundary conditions has the C^^ symmetry, we concen- 
trate on the lowest state with the d^-2_y2 symmetry as 
the two-hole GS. Although this state is not the GS for 
small J (J < 0.43) and Ao = due to finite-size effects, 
this choice is justified by the fact that the two-hole GS 
has dj.2_y2 symmetry in the thermodynamic limit. 
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FIG. 2. Absolute value of the kinetic energy E^in per 
hole for n = 0.2. Open (Solid) symbols correspond to the 
breathing (buckling) mode. Solid (dashed) lines refer to the 
two (one) hole GS. 

For the purpose of our discussion it is convenient to re- 
write the electron-phonon interaction in the boson rep- 
resentation of the phonons, 

E,.pn = ^Y.^h\,h,s^\) (2) 
+Ao5](6i,^+6j,5)(nf T<5) 

i,<5 



where \q = g ■ Since the phononic Hilbcrt space 

has an infinite dimension, we truncate it to a finite num- 
ber of bosonic states i.e. b\^h\^5 < Up^ at each oxygen 
site. We restrict ourselves to Uph = 1. We have tested 
the validity of the one-phonon approximation on the 2x2 
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FIG. 3. Spin structure factor Ssin, n) for a single hole and 
fl — 0.2. Open (Solid) symbols correspond to the breathing 
(buckling) mode. 

As a preliminary study, let us first briefly investigate 
the behavior of a single hole. The absolute values of 
the kinetic energy (t term in Eq. (^) per hole E^in are 
shown as a function of Ao in Fig. ^. For the one- hole GS, 
Ekin decreases significantly with increasing Ao for the 
breathing mode, while it does not change significantly 
for the buckling mode. This is a signature that only the 
breathing mode leads to a polaronic self-trapping process. 
The difference between the two modes is also clear from 
the behavior of the spin structure factor in the one-hole 
GS 

5.(^,7r) = ((^(-l)('==+^»)5f)2) (3) 
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shown in Fig. ||. A significant increase of S's (7r,7r) oc- 
curs around Ao = 0.1 almost independently of J for the 
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breathing mode. The agreement between the behaviors 
of Ss(TT,Tr) and Ekin vs Aq suggests that the increase of 
the effective mass of the hole due to a polaronic self- 
locahzation effect leads, for the breathing mode, to an 
enhancement of the antifcrromagnetic spin correlation. 
Fig. H also shows that the buckling mode, on the con- 
trary, does not lead to any crossover characteristic of self- 
localization. 
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FIG. 4. Two-hol^ binding energy A2 for'' Q, = 0.2 (a) 
plotted vs Ao for fixed J, (b) plotted vs J for fixed A. Open 
(Solid) symbols correspond to the breathing (buckling) mode. 



The two-hole binding energy is a good probe to test 
the formation of pair of holes. It is defined as A2 — 



(2) 



E, 



(0) 



2^;^^^ where E^^"' is the GS energy for a 
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system with Nh = J2i — ^0°'' corresponds to the 
energy of the antifcrromagnetic background. A negative 
value of A2 indicates the stability of a two-hole bound 
state. Fig. ^(a), where A2 is displayed as a function 
of Aq, clearly shows that the buckling mode stabilizes 
the two-hole bound state while the breathing mode sup- 
presses it. The effect of the electron-phonon interaction 
is to shift the boundary of the pairing phase of the t- 
J model: the buckling mode enlarges the phase toward 
small J while the breathing mode reduces it, as revealed 
in Fig. ^(b). The behavior of A2 suggests the possibility 
that the buckling mode assists superconductivity in the 
high- Tc cuprates, while the breathing mode suppresses 
it. We note that, for the buckling mode, no self-trapping 
process occurs even in the two-hole state, since there is 
no significant decrease of the kinetic energy in Fig. |^. 
Thus the hole pair is not localized and can contribute to 
superconductivity. 

The previous data suggest that the electron-phonon 
interaction acts as an effective attraction (repulsion) be- 
tween holes for the buckling (breathing) mode, appar- 
ently in agreement with the anti-adiabatic limit. How- 
ever, for a finite phonon frequency, it is not clear, a priori. 



whether the phonon mediated interaction can be reduced 
to a static potential. To test this possibility, we consider 
the expectation value of the hole-hole distance in the two- 
hole GS 



(4) 



Fig. 1^ shows that, for the breathing mode, dh in- 
creases monotonously with increasing Aq in agreement 
with the effective NN hole-hole repulsion derived in the 
anti-adiabatic limit. However, dh for the buckling mode 
does not show the behavior expected for a NN static at- 
traction. On the contrary, it would rather correspond to 
a small NN static repulsion, at least for small Ao(< 0.2). 
The failure of the anti-adiabatic picture, in this case, sug- 
gests that the effective hole-hole interaction stabilizing 
the hole pairing is controled by an essentially dynam- 
ical effect of the electron-phonon interaction and some 
retardation makes the range of the effective hole-hole at- 
traction longer. In other words, the stabilization of the 
hole binding cannot be understood simply in terms of 
a static NN attraction but rather involves more subtle 
retardation effects. 
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FIG. 5. Hole-hole distance in the two hole GS for S7 — 0.2. 
Open (Solid) symbols correspond to the breathing (buckling) 
mode. 

The dynamical effect of the electron-phonon interac- 
tion can be estimated by the lattice deformation (per 
hole) around the hole sites 



(5) 



where - (-I-) corresponds to the breathing (buckling) 
mode. It is proportional to the absolute value of the 
energy of the electron-phonon interaction. Dh is always 
positive, which means that the oxygen ion deformation 
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toward the neighboring hole sites is favored. Dh is essen- 
tially a dynamical quantity which should be distinguished 
from the total lattice deformaton Dtot given by replac- 
ing {n^ =F n^^g)/Nh by unity in the form (|^). The static 
quantity Dtot is always zero except where $7 = 0. Fig. ^ 
shows that, for the buckling mode, the deformation per 
hole Dh in the two-hole state is larger than the one in the 
single hole state, in contrast to the case of the breathing 
mode. Thus, with buckling modes, the pair takes advan- 
tage of a larger deformation around the holes. On the 
contrary, the breathing deformations lead to an energy 
loss in the pairing state. Thus, the relative change of the 
lattice contraction around the holes in the paired state 
ultimately contribute to a decrease or an increase of the 
binding energy. Particularly for the buckling mode, the 
energy gain coming from the lattice deformation clearly 
dominates the behavior of the binding energy, since no 
significant change appears in the kinetic energy (Fig. |^) 
or in the antiferromagnetic correlation which can be es- 
timated from the spin structure factor 5*8 (tt, tt) in Fig. 
In other words, buckling modes stabilize the hole pair- 
ing state dynamically, with little changes in the static 
features. 



In summary, exact diagonalization studies of the gener- 
alized 2D J-Holstein model give evidence for a stabiliza- 
tion of the two-hole pairing by out-of-plane buckling vi- 
brations of the in-plane oxygens in the high- Tc cuprates. 
On the contrary, in-plane breathing modes suppress the 
pairing. The difference comes from the dynamical effect 
of the lattice displacements which cannot be reduced to 
a simple NN static interaction. In addition, we found 
that the buckling mode does not give rise to any signif- 
icant polaronic self- localization effect, in contrast to the 
breathing mode. 
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FIG. 6. Lattice deformation (per hole) around the hole 
sites for f2 = 0.2. Open (Solid) symbols correspond to the 
breathing (buckling) mode. Solid (dashed) lines refer to the 
two (one) hole GS. 
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This is to be contrasted to the case of the breathing 
mode which can modify some static properties. As Aq in- 
creases, the effective mass becomes larger due to the po- 
laronic self-trapping process and, as a result, the antifero- 
magnetic correlation increases. Since the self-localization 
effect in the two-hole state is smaller than the one-hole 
state (Fig. ||), this effect might tend to stabihze some 
trapped two-hole bound states. However, a larger effec- 
tive repulsion due to the dynamical lattice deformation 
overcomes the static effect and suppresses the hole pair- 
ing in the presence of breathing modes. 
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